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We present our most recent results on the charmonium spectrum using relativistic Wilson fermions. We study 
the dependence of the spectrum on the charm quark mass and the Wohlert-Sheikholeslami improvement term. 
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1. INTRODUCTION 

This is a status report on our calculation of 
the charmonium spectrum. Our first results were 
presented last year ||l|. An important aspect of 
this calculation, the determination of as from the 
charmonium spectrum, has already been reported 
separately The detailed results of this study 
will be presented in a future publication ^ . 

The calculation of experimentally well mea- 
sured quantities in the charmonium spectrum 
serves several purposes. The investigation of 
physical quantities that are sensitive to the sys- 
tematic errors present in lattice calculations can 
lead to a better understanding and possibly the 
removal of these errors. On the other hand, a 
quantity that is relatively independent of the sys- 
tematic errors of the lattice calculation can be 
used to obtain information about fundamental 
parameters of the Standard Model. 

In Refs. ||l|,|| we argued that the IP- IS splitting 
falls into the latter class of quantities and there- 
fore lends itself to an accurate determination of 
the strong coupling constant of QCD. The split- 
ting between the J/^ and 77c, on the other hand, 
is expected to strongly depend on the leading lat- 
tice spacing errors. This splitting and the Icptonic 
decay amplitude of the J/ip are expected to also 
show some dependence on the charm quark mass. 
One might therefore use these two quantities in a 
study of the former type. 

Some of the details of the calculation are de- 
scribed in section ||. The effect of the omission of 
sea quarks has already been discussed in detail in 
our analysis of as ■ In section ^ we discuss their 
effect on the hyperfine splitting and the leptonic 



decay amplitude of the J/ 4'- The results of our 
study of the dependence of the spectrum on the 
charm quark mass and 0{a) lattice spacing errors 
are presented in section ^. 

2. THE CALCULATION 

We analyze three different lattices (12'^ x 24, 
16^ X 32, 24"*) at three different couphngs (/3 = 
5.7,5.9,6.1) such that the volumes are similar 
while the lattice spacing varies by about a factor 
of two (see Ref. ||l|,^ for more details). In order 
to investigate the effects of lattice spacing errors 
on the spectrum, we use Wilson fermions with 
(c — 1.4) and without (c = 0) the 0{a) correction 
term . We study the charm quark mass depen- 
dence by varying the charm hopping parameter, 
Kcharm, over a reasonable range for each lattice 
and fcrmion action (or choice of c). 

The 2-pt. functions used for the mass splittings 
have already been described in detail in Ref. [Q. 
The mesons are created and destroyed by oper- 
ators, x(a;) = X)r''/'(-^ + '')r'0(a;)S'(r), with expo- 
nential spreading functions, S'(r), that optimize 
the overlap with the ground states. The leptonic 
matrix element for the J/ip is extracted from the 
2-pt. function (F = 7^) 
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V^{x) — %l!{x)^^ij}{x) is the local vector current. 
With the convention {p\p') = {2ttY6'^{p — p'), 
Eq. (P takes the asymptotic form 



\V^\J/i,){J/i,\x^\Q)e' 



(2) 
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We parametrize the matrix element in terms of 
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{Q\V^\J/^)=e^V^ 



(3) 



is related to the decay constant ^ used by 

other groups 1 by = (m^ /2)i/2/^'. We 
renormalize the current non-perturbatively with 
the charge matrix element {J / ip\V4\J / ^) = Zy^. 



As remarked in Refs. ||6|,[7|, the perturbative cal- 
culation of the local current renormalization Zy 
(using improved perturbation theory repro- 
duces the non-perturbative result to a few % in 
the case of Wilson (c = 0) fermions, if the quark 
fields are properly normalized. 

3. THE EFFECT OF SEA QUARKS 

A nonrelativistic system, like charmonium, is 
well described by a potential. We can therefore 
understand the dominant effects of the sea quarks 
on the spectrum via their effect on the poten- 
tial. Following the arguments of Ref. in set- 
ting the scale with the IP-IS splitting, we 
require in effect the potentials of the full and the 
quenched theory to match at the relevant physics 
scale which is the intermediate distance scale of 

400 - 750 MeV. We then expect the (perturba- 
tive) short distance potential to be too weak, as a 
result of the zero flavor /3-function being slightly 
too large. 

For quantities that are dominated by short dis- 
tance effects, like the hyperfine splitting and lep- 
tonic decay amplitude of the J/tp: we have to con- 
sider the sea quark effects on the short distance 
potential. This is best investigated within the 
context of a potential model. One expects for the 
hyperfine splitting to lowest order 0] 



Am(J/V' — rjc 
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*(0)p 



(4) 



Similarly, according to the van Royen - Weisskopf 
formula we expect for the leptonic matrix element 



l*(0)| 



(5) 



Thus, the effect of the sea quarks on the hyperfine 
splitting and the leptonic matrix element can be 
estimated by their effect on the wave function at 
the origin, ^'(0), in a potential model. 



The Richardson potential is convenient for 
our purposes, because it incorporates asymptotic 
freedom and confinement in a simple form [|fO|: 



with /3^"^^ = ff - 2?i//3. Figure g shows the 
wave functions of the IS state, obtained from fit- 
ting the Richardson potential to the experimental 
charmonium spectrum with Uf = 3 and rif = in 
comparison. We find for the ratio of wave func- 
tions at the origin: 
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g2 ln(f + 9VA2) 



(6) 



^.(o)(0) 



= 0.86 



(7) 



^'(3)(0) 

The hyperfine splitting receives an additional 
correction from as{mc). The calculation of the 
effects of the sea quarks on as{7r/a) is described 
in detail in Ref. The resulting correction for 
as{mc) is analogously: 

(0)/ ^ 



a. 



(3), 



= 0.8f ±0.06 



(8) 




Figure f. The radial wave function for the fS 
state in the Richardson potential model. The 
solid line is for Uf = 3 and the dashed line for 
n/ = 0. 
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In summary, we estimate the quenched hyperfine 
spHtting to be reduced by 40 % from the experi- 
mental value, Ato(J/-0 - 7]^^^ = 117.3 MeV to 



Am( J/V' - ry,)q"™=hod _ 



(9) 



The quenched leptonic matrix element just re- 
cieves the correction in Eq. (|^), a 14% reduc- 
tion from the experimental result V^,^^ = 0.509 
GeV3/2: 



^quenched ^ Q^^gg GeV^^^ 



(10) 



There is certainly no reason to believe that 
the sea quarks do not change the potential (non- 
perturbatively) in a way that has not been taken 
into account here. Eqs. and ( |lO| ) are there- 
fore only rough estimates of the possible effect. 
However, from the above considerations it should 
be expected that the omission of sea quarks re- 
duces both the hyperfine splitting and the lep- 
tonic width from their experimental values. 

4. RESULTS 

Figures |, | and | show the IP-IS sphtting, 
the J/ijj - rfc splitting, and the leptonic matrix el- 
ement of the J I i/j as functions of the charm hop- 
ping parameter for Wilson fermions with (c — 
1.4) and without (c — 0) the improvement term 
on the 16'^ x 32 (/3 = 5.9) lattice as representative 
examples of our results on all three lattices. As 
expected, we find the dependence of the IP-IS 
splitting on the charm quark mass (parametrized 
by Kcharm) and the coefficient of the improvement 
term, c, to be very weak, smaller than our statisti- 
cal errrors. This confirms our previous arguments 
|lj0] to use this splitting to determine the scale 
a in lattice calculations and subsequently ex- 
tract the strong coupling constant from it. Our 
previous result for remains unchanged after 
this study. 

Both the hyperfine splitting as well as the lep- 
tonic matrix element vary significantly with the 
two action parameters, Kcharm and c. However, 
from the considerations in section |[ it is clear 
that a phenomenological determination of these 
parameters is not possible within the quenched 
approximation without a better understanding of 




Figure 2. The IP-IS splitting in lattice units vs. 
Kcharm On the 16"^ X 32, /3 = 5.9 lattice; the circles 
arc for c = 1.4, the squares are for c — 0. 
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Figure 3. The J/tp-ijc splitting in lattice units vs. 
«^charm On the 16"^ X 32, P = 5.9 lattice; the circles 
are for c — 1.4, the squares are for c = 0. The 
dashed line is the experimnetal splitting; the solid 
line is the splitting of Eq. (9) . 
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Figure 4. The leptonic matrix element, K/i in lat- 
tice units, vs. Kcharm On the 16*^ x 32, (3 — 5.9 
lattice; the circles are for c = 1.4, the squares are 
for c = 0. The dashed line is the experimental 
value; the solid line is the value of Eq. (10). 
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the effects of the sea quarks on the static po- 
tential. We therefore take the mean field value 
c = 1.4 as our best estimate for the coefficient of 
the improvement term. Setting the charm quark 
mass with the leptonic matrix element, (us- 



ing Eq. (10)), we extract hyperfine splittings as 
shown in figure || for all three lattices. The split- 
ting in physical units is graphed as a function of 
for better visiblility. Also shown are the exper- 
imental and the expected quenched values. We 
find very little variation of the hyperfine splitting 
with the lattice spacing, indicating small residual 
lattice spacing errors. Our final result 



Am( J/^ — rjc 



93 ± 10 MeV 



(11) 



(statistical error only) lies within the expectations 
for a quenched calculation, slightly below the ex- 
perimentally measured splitting. 
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Figure 5. The J/tp-rjc splitting vs. a^. The two 
points on the left (burst) represent the experi- 
mental splitting (upper) and the expectation for 
the quenched result (lower). 
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